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Abstract 

We have measured the collective excitation spectrum of interacting electrons in 
one-dimension. The experiment consists of controlling the energy and momentum 
of electrons tunneling between two clean and closely situated, parallel quantum wires 
in a GaAs/AlGaAs heterostructure while measuring the resulting conductance. We 
measure excitation spectra that clearly deviate from the non-interacting spectrum, 
attesting to the importance of Coulomb interactions. Notable is an observed 30% 
enhancement of the velocity of the main excitation branch relative to non-interacting 
electrons with the same density. In short wires, finite size effects resulting from 
broken translational invariance are observed. Spin - charge separation is manifested 
through moire patterns, reflecting different spin and charge excitation velocities. 
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The problem of interacting electrons in one dimension has attracted consid- 
erable attention over the years, both theoretically and experimentally. One 
dimensional systems differ fundamentally from their higher dimensional coun- 
terparts in the role interactions play. The latter are well described by Landau 
Fermi-liquid theory [1], in which a system of interacting electrons is mapped 
onto a system of weakly interacting long-lived quasiparticles [2] , each carrying 
a fundamental unit of charge, — e, and a fundamental unit of spin, h/2. In 
one dimension the quasiparticles become unstable because of their large prob- 
ability to decay into particle-hole pairs, a process that is blocked in higher 
dimensions due to phase-space constraints. A ID system is best described as 
a Luttinger-liquid [3,4], characterized by decoupled long-range spin and den- 
sity correlations. It is predicted that an injected electron will break up into 
separate spin and charge excitations which propagate with separate velocities 
and give separate singularities in the one-electron Green function. 

In the experimental study of clean one-dimensional (ID) systems we search 
for effects that are a direct and exclusive consequence of interactions and the 
correlations they induce. The first transport studies of wires aimed at measur- 
ing the conductance G = dl/dVso (I is the current along a one-dimensional 
system subjected to a bias V S d) [5,6]. Only later was it realized that in clean, 
finite wires G is quantized in units of e 2 /h regardless of the strength of in- 
teractions [7,8]. The next step was the realization that disorder allows to tap 
into the unique properties of a Luttinger-liquid. In the limit of low energies, 
a single, weak scatterer in a ID wire with repulsive interactions effectively 
cuts it into two disjoint pieces. Thus G reflects the long range correlations in 
a Luttinger-liquid through a power-law dependence on energy [9,10,11,12,13], 
that is absent in a Fermi-liquid. 

When looking for ways to study a many body system the most fundamental 
quantity to focus on is its excitation spectrum. Naturally one would espe- 
cially like to be able to measure it directly. First steps in this direction were 
experiments that used either Raman spectroscopy [14] or angle resolved pho- 
toemission [15,16,17,18] to measure the excitation spectrum of electrons in 
one dimension. Below we show that measuring tunneling between two par- 
allel wires enables to circumvent some of the limitations inherent in optical 
spectroscopy [19,20] and allows to map out the dispersions of electrons in one 
dimension up to and beyond the Fermi momentum with high accuracy. 

The main idea behind using tunneling amongst parallel wires for spectroscopy 
is that it affords simultaneous control of both the energy and the momentum 
transferred between the wires by the tunneling electrons. The energy is given 
by eV, where V is the bias voltage difference between the wires while the 
momentum is hqs = eBd, where B is a magnetic field applied perpendicular 
to the plane of the wires and d is the distance between them. For long junc- 
tions the tunneling process conserves both energy and momentum, so unless 
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Fig. 1. Schematic of the sample and the contacting scheme. The sample is fabricated 
using CEO. Two parallel ID wires (dark gray) span along the whole cleaved edge. 
The upper wire overlaps the 2DEG (light gray) in the upper quantum well (20nm 
thick), while the lower wire is separated from them by a 6nm insulating barrier in 
an otherwise empty quantum well (30nm thick). Contacts to the wires are made 
through the 2DEG. Several 2//m-wide tungsten top gates can be biased to deplete 
the electrons under them (only G\ and G2 are shown). The magnetic field B is 
perpendicular to the plane defined by the wires. The depicted configuration allows to 
study the conductance of a single wire- wire tunnel junction of length L by measuring 
the current / that flows when a bias voltage V is applied between the wires. 



there exists a state in the target wire at the particular value of energy and 
momentum selected by V and B, tunneling is blocked. This allows to directly 
determine the dispersions of the many-body states (see e.g. [21,22]). 

The first study of tunneling between two parallel wires was facilitated by 
cleaved-edge overgrowth [23,6,24] of a AlGaAs/GaAs double quantum well 
heterostructure (see Fig. 1). An initial growth sequence renders the 20nm- 
wide upper quantum well occupied by a two-dimensional electron gas (2DEG) 
and the 30nm-wide lower quantum well devoid of electrons. The two quantum 
wells are separated by a 6nm-wide insulating AlGaAs barrier. After a 30sec 
infra-red illumination the mobility of the 2DEG is fi ~ 3 x 10 6 cm 2 V s _1 and 
its density is n ~ 2 x 10 11 cm -2 . After depositing 2/im tungsten gates on the top 
surface of the sample, it is reinserted into the molecular-beam epitaxy (MBE) 
chamber. The sample is then cleaved in this pristine environment to expose 
a clean (110) plane and a second growth sequence is initiated. As a result a 
trapping potential for electrons is created along the cleaved edge in each of the 
wells. The quantum states in this potential constitute the sub-bands in each 
of the wires in the experiment. Typically there are 3-5 occupied sub-bands in 
each of the wires following the 30sec illumination mentioned above. 

The measurements reported here were conducted in a 3 He refrigerator at a 
base temperature of 0.25K (unless otherwise stated) using a lockin amplifier 
at a frequency of 14Hz with an excitation of 10/xV. All of the measurements 
were two-terminal measurements between indium contacts to the 2DEG in the 
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Fig. 2. Plot of G (V, B) for a 10//m junction. In order to highlight the features a 
smoothed background has been subtracted from the raw data. The dispersions of 
the upper wire are easily discernable. Also plotted are fits to the free electron model 
with m* = 0.75m GaAs for the lower wire and m* = 0.85m GaAs for the upper wire 
(corresponding to g v = 0.85 and g L = 0.75). The extracted (zero field) densities are: 
{n L } = 88, 47, 32, 24, 15 ± l^m^ 1 for the lower wire and {%} = 90, 64, 62 ± l/mi^ 1 
for the upper wire. 

upper quantum well (see Fig. 1). The source is the 2DEG between gates G\ 
and G 2 . The bias on G\ is set to deplete both wires and the 2DEG, while the 
bias on G 2 is set to leave only the lower wire conducting. Thus, the upper wire 
between G\ and G 2 is at electrochemical equilibrium with the source 2DEG, 
while the whole semi-infinite lower wire is in equilibrium with the drain, the 
2DEG beyond G 2 . Thus, any voltage difference induced between the source 
and the drain drops on the tunnel junction between G\ and G 2 . 

A typical measurement of G (V, B) for a long junction is shown in Fig. 2. 
The most prominent features in the figure are parabolic-like curves. In the 
noninteracting electron picture, these curves would mark boundaries in the 
V — B plane between regions where tunneling is blocked and regions where it 
is allowed. As explained in [20], each of these curves is the dispersion of a mode 
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in one of the wires. With the voltage convention in Fig. 1, curves with turning 
points at V > are lower wire dispersions and those with turning points 
at V < are upper wire dispersions. At V = only electrons at the Fermi 
level can participate in tunneling. Therefore tunneling is appreciable only if the 
momentum transferred to a tunneling electron is given by ± (kp ± kj^j , where 

kp' L = nn U]L /2 in a spin-degenerate mode (n UjL denote electron densities of 
modes in the upper and lower wires and Zeeman splitting is ignored because it 
is not important here). This happens only when the magnitude of the magnetic 
field is: 



ed 



k u F ±k L F . (1) 



Thus B are a direct measure of the sum and difference of the densities of 
each of the modes in each of the wires. This allows to calculate the dispersions 
non-interacting electrons would have at the same density, because they depend 
only on the band mass m GaAs and the parameters of the heterostructure, which 
are known. We find that the observed dispersions deviate significantly from 
the non-interacting dispersions. 

For quantitative comparison with non-interacting dispersions one has to ac- 
count for the fact that in finite width wires, like those reported here, the 
dispersion of each mode depends on magnetic field. As B is increased the 
dispersions become more flat (when B is extremely large Landau levels are 
be recovered). This affects the densities, the general trend being that lower- 
energy modes are populated at the expense of higher-energy modes. Solving 
the Schrodinger equation numerically for a finite square potential well in the 
growth direction, we determine the zero field occupations from the B ±, s and 
find the shape of the curves traced out in the V — B plane by the dispersions. 
When this is done we find a systematic discrepancy - the real traces have an 
enhanced curvature. To model this observation we solve the finite well problem 
again, but with a renormalized mass, m* = gm GaAs , where < g < 1. This 
corresponds to a dispersion with the same density as for m GaAs , that has an 
enhanced curvature and hence an enhanced Fermi- velocity, vp/g. In all of the 
cases we have studied, g is significantly suppressed from 1. For example, in 
Fig. 2 we find for the lower wire g v = 0.75 and for the upper wire g v = 0.85, 
whereas in Fig. 3 we find g L = g v = 0.7. While the understanding of the cur- 
vature enhancement is still lacking, the enhanced velocity is consistent with 
Luttinger-liquid theory. The theory predicts low voltage peaks in the tunneling 
conductance that trace out dispersions of charge modes that have a velocity 
v c = v F /g [21]. 

The high resolution scan of a 2/zm junction shown in Fig. 3 reveals inter- 
esting extra features. Zooming on the low field region the peaks tracing the 
dispersions can be seen clearly. We have overlayed the plot with our simple cal- 
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Fig. 3. High resolution plot of G (V, B) for a 2//m junction. In order accentuate 
weak features a smoothed background has been subtracted, and the color-scale 
made nonlinear. Light shows positive and dark negative signal. The black curves 
are the expected dispersions of noninteracting electrons at the same electron density 
as the lowest-energy ID bands of the wires, as determined from the crossing points 
with the 5-axis (only the lower one is shown in the figure). The white curves are 
generated in a similar way but with a renormalized GaAs band-structure mass: 
m* = 0.7m GaAs . This corresponds to g v = g L = 0.7. Only the curves labelled a,b,c 
& d in the plot are found to trace out experimentally-observable peaks in G (V, B) 
with the curve d following the measured peak only at V > — lOmV. 

culation of the dispersions: In black we plot the dispersions with m* = m GaAs 
while in white we plot the dispersions obtained with m* = 0.7m GaAs . This 
value of m* was chosen to fit the observed dispersions also near B + (not 
shown in the figure). Lines a,b & d describe the dispersions rather well, at 
least above — lOmV, leading us to conclude that both wires possess modes 
with v c = vp/0.7. Interestingly there is an extra mode in the upper wire, 
marked c that lies very close to a noninteracting curve belonging to the upper 
wire. We can thus conclude that there is a mode in the upper wire that moves 
at approximately vf, as expected in Luttinger-liquid theory for a spin mode. 

In addition to the dispersions seen in Fig. 3 this figure contains features that 
we attribute to the finite length of the junction. The first of these is an intricate 
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pattern of oscillations, to which we return below. The second is a zero bias 
anomaly that is seen as a sharp crease near V = 0. Fig. 4 shows that this dip 
is very sensitive to temperature. 




0.24 0.4 0.8 1.6 

T[K] 



Fig. 4. Conductance near V = as a function of temperature at B = 2.5T (•) for 
a 6/xm junction. The data was extracted from scans such as the those shown in the 
insets. Also shown are the fits to T a ™ d ^ (solid line) and T ab "^ (dashed line), 
where we used g = 0.59. Insets: Non-linear tunneling conductance at B = 2.5T as 
a function of V for T = 0.24K and T = 0.54K. Here we show fits to Eq. 2 with 
a(T) = a cnd (solid line) and a(T) = a bulk (dashed line). 



The source of the dip is the suppression of the tunneling density of states 
characteristic of a Luttinger liquid [25]. For the range of B in Fig. 3, the signal 
is dominated by tunneling from the Fermi-liquid 2DEG (having a constant 
density of states at the Fermi level) in the upper well to a Luttinger-liquid in 
the lower well. The zero-bias dip is thus entirely due to electron correlations 
in the lower wire. This is manifested by [26]: 

G(KT)«r%(^), (2) 

where F a is a known scaling function obeying F a (x) ~ 1 as i -> and 
F a (x) ~ x a for x ^> 1. The exponent a obeys [25]: 
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a(T) ~ < 



"bulk k B T > hv F / (gL) 
ttcnd k B T <^hv F /{gL). 



ttbuik.cnd are the exponents obtained for tunneling into the middle and into the 
end of a Luttinger liquid [27]: 



a bulk = (g + g- 1 -2)/4, (4) 
a cnd = (^ 1 -l)/2. (5) 

In [25] it is argued that a bulk occurs because when T is not too small (relative to 
VF/(gL)), the tunneling process is insensitive to its occurring near the end of 
the lower wire, so the exponent a is characteristic of tunneling into the middle 
of a Luttinger liquid (4). For lower values of T the tunneling is effectively into 
the end of the lower wire so the exponent is the exponent for tunneling into 
the end of a Luttinger liquid (5). 

To compare the data in Fig. 4 to Eq. 2 we first fit the data to the V = limit 
of Eq. 2: 

G ~ T a<yT \ (6) 



where ot(T) is given by Eq. 3. The result of the fit is overlayed on the data in 
Fig. 4. As a corroboration we use the parameters from the fit in Eq. 2 with 
either a bulk or a cnd . The results are plotted in the insets of Fig. 4, where one can 
see that at low T and V the data is reasonably well described by Eq. 2 with 
ot(T) = a cnd while when T or V increase there is a crossover and a{T) = a bulk . 

We now return to the oscillation pattern seen in Fig. 3 and in Figs. 5, 6. In 
the last two figures the range of field is such that the lines that correspond 
to the dispersion curves appear as pronounced peaks that extend diagonally 
across the figures. In addition to these we observe numerous secondary peaks 
running parallel to the dispersions of the lowest-energy modes. These side 
lobes always appear to the right of upper wire dispersions, in the region that 
corresponds to momentum conserving tunneling for an upper wire with a 
reduced density. Thus, the V — B plane is separated into quadrants: quadrant 
I has a checkerboard pattern of oscillations, quadrant II has a hatched pattern 
and quadrant III has no regular pattern. Furthermore, the frequency of the 
oscillations depends on the length of the junction. When the lithographic 
length of the junction, L, is increased from 2/im (Fig. 5) to 6/im (Fig. 6) the 
frequency in bias and in field increases by a factor of approximately three. 
The period is related to the length of the junction, L, by the formula: 

AVL/v p = ABLd « 0o, (7) 
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Fig. 5. (a) Conductance oscillations at low field from a 2//m junction. A smoothed 
background has been removed from the raw data and the nonlinear scale has been 
optimized to increase the visibility of the oscillations. Parallel side lobes attributed 
to finite-size effects appear only to the right of the main dispersion peaks, defining 
quadrants I, II and III. Also present is a slow modulation of the interference along 
the F-axis. (b) Absolute value and position (converted to length) of the peak cor- 
responding to the oscillations along B in quadrant II of (a), as determined from the 
Fourier transform of S^^G (V, S^ 1 ^). See main text for definition of S, (3 and 
other details. The slow modulation as a function of V is easily discerned. 

where 0o = h/e is the quantum of flux and v p is an effective velocity. Eq. 7 
can be used to extract L, because d is known. 

To understand the asymmetry of the interference around the dispersion curves 
it suffices to consider non-interacting electrons in both wires. Then, for weak 
tunneling, the current is related by the Fermi golden rule to the tunneling 
matrix element between a state in one wire and a state in the other [28,25]: 

I(V,B) ~V \\M{k + )\ 2 + \M{kJ)\ 2 ] , (8) 
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where 
M 



OO 

k) = J dxe iKX i)jj(x)e~ ik F x (9) 



and k± = kp — kp + eV/(hvp) ± Here vf is the Fermi velocity, which 
is almost the same in both modes giving rise to the interference. We shall 
further assume that the potential limiting the upper wire's length, U(x), is 
smooth on the scale of the Fermi wavelength, a reasonable assumption since 
it is defined by top gates lying 0.5/im away. Under this assumption the WKB 
approximation can be used to write: 

ipu{%) ~ k~ l l 2 (x) exp ikpX — is(x) , (10) 



where k(x) = kpJl - U(x)/E^ and s(x) = Jq dx' kp - k(x' 



Eq. 8 can be 



used to find U(x). We have found that to a good approximation it is given 
by U(x) « Ep \2x/L eS f. In this model, the function S^^G (V, S 1+1 ^) is 

periodic in S' 1+1///3 , with a period determined by L eB , and where S = hn + /(ed) 
is in quadrant II in Figs. 5a and 6a. We find that the data is well described by 
f3 — 8 ± 2 for the short junction and (3 = 21.5 ± 2 for the long junction [25]. 
With these values of f3 we then perform Fourier analysis of the data for each 
value of V . The results for the data in Fig. 5a are plotted in Fig. 5b, were 
both the amplitude of the main peak and its position are shown. In the figure 
we convert the frequency of the oscillations back to the effective length for 
the upper wire, L eS (V}, which clearly depends only very weakly on V. Similar 
analysis was performed for the data in Fig. 6a. The length extracted from the 
Fourier analysis for both the 2/xm junction (L cS — 2.7 ± 0.1/xm) and for the 
6/xm junction (L eff = 7.3 ± 0.3/xm) is larger than the lithographic length. This 
is reasonable because the gates delimiting the upper wire are on the surface 
of the sample. 

Next we turn our attention to the structure of the amplitude of the main 
Fourier peak (see panels b in Figs. 5, 6). The amplitude is seen to oscillate as a 
function of V, giving rise to a series of vertical strips of suppressed conductance 
in Figs. 5a and 6a. To understand this effect the finite interactions in the wires 
have to be taken into account. A more general expression for the tunneling 
current (to lowest order in tunneling) is given by [21]: 

oo oo oo 

I{V,B)oc J dx J dx' J dte iqB{x - x,) e ieVt/n C(x,x';t). (11) 



where C(x, x';t) is the two point Green function, which is known from Lut- 
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Fig. 6. Same as Fig. 5 but for a junction. Note that the oscillations are approx- 
imately three times faster than in Fig. 5, as expected from Eq. 7. For this junction, 
the asymmetry in the strength of the side lobes on opposite sides of a dispersion 
peak is less pronounced than for the shorter junction appearing in Fig. 5. 



tinger liquid theory. In the case we are describing here one finds that the 
interference has a contribution from two velocities in the upper wire, vp and a 
charge mode velocity, v~. v~ is the velocity of the antisymmetric charge mode 
which arises because the Fermi velocities in the two wires are very similar. 
As a result of this, the interference in Figs. 5 and 6 can be understood as a 
moire pattern created by side lobes running with a slope (in the V — B plane) 
{ypd)~ l and other side lobes with a slope (t>~<f) _1 . From the ratio between AV 
(denned in Eq. 7) and the distance between the suppression strips, AV mod , the 
ratio between the two velocities can be extracted. One finds: 



AV 
AV 



11 + 9- 



21 



9- 



(12) 
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where g_ = vf/v ~. For the data presented here we find g_ = 0.67±0.07. This 
is in agreement with previous assessments of g in CEO wires and is a direct 
consequence of spin-charge separation in our wires. 
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